gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十六章/eg16_2/eg16_2_1.m

    x=1:0.1:20;
disN=40/pi./x.*sqrt(1-x.^(-2))-j*40/pi./x.^2;                                   %描述函数
disN2=-1./disN;                                                                            %负倒描述函数
w=1:0.01:200;
num=12;                                                                                       %线性部分分子
den=conv([1 1],[1 6 13]);                                                              %线性部分分母
[rem,img,w]=Nyquist(num,den,w);                                               %线性部分Nyquist曲线参数
plot(real(disN2),imag(disN2),rem,img)                                         %同时绘制非线性部分和线性部分的极坐标图
grid; xlabel('Re');ylabel('Im');
%读出线性部分和非线性部分交点的坐标值，并利用坐标值求出振荡幅值和频率
w0=spline(img,w,-0.0785)                               %当img=-0.0785时，所对应的w值
x0=spline(real(disN2),x,-0.166)                       %当disN2的实部为-0.166时，所对应的x值